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Abstract 

Specialized microenvironments called niches regulate tissue homeostasis by controlling the balance between stem cell self- 
renewal and the differentiation of stem cell daughters. However the mechanisms that govern the formation, size and 
signaling of /n vivo niches remain poorly understood. Loss of the highly conserved histone demethylase Lsdl in Drosophila 
escort cells results in increased BMP signaling outside the cap cell niche and an expanded germline stem cell (GSC) 
phenotype. Here we present evidence that loss of Lsdl also results in gradual changes in escort cell morphology and their 
eventual death. To better characterize the function of Lsdl in different cell populations within the ovary, we performed 
Chromatin immunoprecipitation coupled with massive parallel sequencing (ChlP-seq). This analysis shows that Lsdl 
associates with a surprisingly limited number of sites in escort cells and fewer, and often, different sites in cap cells. These 
findings indicate that Lsdl exhibits highly selective binding that depends greatly on specific cellular contexts. Lsdl does not 
directly target the dpp locus in escort cells. Instead, Lsdl regulates engrailed expression and disruption of engrailed and its 
putative downstream target hedgehog suppress the Lsdl mutant phenotype. Interestingly, over-expression oi engrailed, but 
not hedgehog, results in an expansion of GSC cells, marked by the expansion of BMP signaling. Knockdown of other 
potential direct Lsdl target genes, not obviously linked to BMP signaling, also partially suppresses the Lsdl mutant 
phenotype. These results suggest that Lsdl restricts the number of GSC-like cells by regulating a diverse group of genes and 
provide further evidence that escort cell function must be carefully controlled during development and adulthood to ensure 
proper germline differentiation. 
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Introduction 

Stem cells undergo self-renewing divisions in which at least one 
daughter retains its stem cell identity, while the second daughter 
may or may not difiFerentiate, depending on intrinsic and extrinsic 
cues. A balance between stem cell self-renewal and difiFerentiation 
must be maintained for proper organ formation during develop- 
ment and tissue homeostasis in adulthood. Stem cells often reside 
in microenvironments called niches, and specific mechanisms 
tightly regulate the size and signaling output of these structures [1]. 
However, in vivo niches have often proven difficult to identify in 
mammalian tissues. As a result, much of the current understanding 
of niches stems from the study of invertebrate models such as the 
germline stem cells (GSCs) of the Drosophila ovary. 

Drosophila female GSCs reside in a well-characterized niche at 
the tip of a structure called a germarium (Figure lA). Within 
germaria, GSCs lie immediately next to a somatic cell niche 
comprised of cap cells and terminal filament cells [2] . Escort cells 
reside adjacent to the cap cells and line the anterior portion of the 
germarium. These cells act to shepherd the germ cells during the 
earliest stages of their differentiation [3,4], after which developing 



germline cysts are enveloped by follicle cells derived from a second 
stem cell population within the germarium. Cap cells produce 
Decapentaplegic (Dpp), which in turn activates a canonical Bone 
Morphogenic Protein (BMP) signal transduction pathway in GSCs 
[5,6]. BMP pathway activation results in the transcriptional 
repression of hag of marhles (ham) [7-9] , a factor both necessary and 
sufficient for germ cell dilferentiation [10,11]. Ectopic Dpp 
signaling outside the tip of the germarium results in an expanded 
GSC phenotype [5,9]. 

Other pathways and neighboring cells likely regulate niche 
specific BMP signaling. For example, a recent study provides 
evidence that hedgehog [hh) produced by the cap cells stimulates the 
anterior escort cells to produce niche specific signals [12] 
Moreover, several additional intrinsic and extrinsic mechanisms 
help restrict Dpp ligand production and diffusion within the niche 
(reviewed in [13,14]). One such mechanism involves the histone 
demethylase Lysine Specific Demethylase 1 (Lsdl). Lsdl uses a 
flavin-dependent monoamine oxidative based mechanism to 
remove mono- and di-methyl groups from histone H3 on lysine 
4 (H3K4mel and H3K4me2) [15]. In mammals, Lsdl has been 
shown to silence a number of distinct gene sets in different cellular 
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Author Summary 

The mechanisms that govern the formation, size and 
signaling output of in vivo niches remain poorly under- 
stood. Studies of Drosophila germline stem cells (GSCs) 
have suggested that chromatin programming greatly 
influences the behavior of these cells and their progeny. 
Previous work has shown that loss of the highly conserved 
histone demethylase Lsdl results in ectopic niche signal- 
ing and an expanded GSC phenotype. To determine direct 
regulatory targets of Lsdl, we employed chromatin 
immunoprecipitatlon coupled with massive parallel se- 
quencing (ChlP-seq) using specific cell populations inside 
and outside of the GSC niche. These experiments revealed 
that Lsdl exhibits highly enriched binding to over one 
hundred genomic sites within a specific cell population. 
Furthermore, mis-regulation of some of these direct 
targets contributes to the expanded stem cell phenotype 
observed in Lsdl mutants. These results provide insights 
into how Lsdl directly restricts the size of the GSC 
microenvironment and establish a platform for under- 
standing and exploring chromatin programming inside 
and outside an in vivo stem cell niche. 



contexts, including Notch targets, TGFP-1 and various loci 
involved in the maintenance of embryonic stem cells [16-20]. 
Additional studies suggest that Lsdl may also promote gene 
expression under certain circumstances [21]. Disruption of 
Drosophila Lsdl results in a male and female sterUe phenotype, 
marked by the expansion of GSC-like cells in the germarium 
[22,23]. These cells exhibit ectopic BMP responsiveness and fail to 
initiate a normal differentiation program once they leave the cap 
cell niche [24] . 

To characterize the molecular mechanism by which Lsdl 
restricts signaling outside the Drosophila female GSC niche, we used 
ChlP-seq to define direct binding sites of Lsdl specifically in either 
escort cells or cap cells. These experiments revealed that Lsdl 
binds to over one hundred sites in escort cells and provide further 
insights into how Lsd 1 contributes to the chromatin programming 
of cells inside and outside of an in vivo niche. 

Results 

Further characterization of the Lsdl loss-of-function 
phenotype 

Escort cells send out extensions that closely contact germline 
cysts undergoing the early steps of differentiation [3,4] . Escort cell 
death or genetically disrupting escort cell extensions can lead to an 
inappropriate expansion of GSC-like cells in the germarium [3] . 
Previous results showed that Lsdl functions within escort cells to 
prevent expanded BMP signaling outside of the GSC niche [24] . 
This phenotype was accompanied by widespread cell death in 
both somatic cells and germ cells. Therefore we considered the 
possibility that the expansion of BMP signaling exhibited by Lsdl 
mutants may depend on changes to escort cell morphology and 
number. To test this, we knocked down the expression of Lsdl 
specifically in the escort cells and early follicle cells by crossing 
UAS-Lsdl^"^' into a c587-gal4; UAS~mCD8::GFP background and 
stained the resulting ovaries for GFP and the fusome marker Hts. 
Fusomes are highly vesiculated organelles that appear round in 
GSCs and cystoblasts, and become branched as these germ cells 
differentiate into multi-cellular cysts [25,26]. Three days after 
eclosion, control samples appeared normal. These germaria 
typically contained two to four single cells (GSCs and cystoblasts) 



with round fusomes and escort cells that extended cytoplasmic 
processes between the developing cysts (Figure IB). In contrast, the 
Lsdl RNAi samples showed an expansion of GSC-like cells with 
round fusomes. Escort cell extensions were clearly present in some 
germaria, but were missing in others (Figure 1C,D). These 
observations suggested that while the knockdown of Lsdl caused 
changes in escort cell morphology, the presence of extra single cells 
with round fusomes did not absolutely depend on a complete loss 
of escort cell extensions. However changes in escort cell 
morphology likely contributed to the phenotype over time. In 
addition, expression of Lsdl^'^' also led to an increase of death 
within escort cells, consistent with the widespread cell death 
previously noted in Lsdl nuU mutant germaria (Figure 1 E,F) [24] . 

Next we performed clonal analysis using the mosaic analysis 
with a repressible cell marker (MARCM) system to further analyze 
the Lsdl mutant phenotype. Clonal germaria were stained for the 
positive clone marker GFP and for the fusome marker Hts. We 
categorized the relative position of escort cell clones along the 
anterior-posterior axis of the germarium. Cells were considered 
anterior escort cells if they were immediately next to the cap cells, 
posterior escort cells if they were immediately adjacent to the 
follicle stem cells and middle escort cells if they were located in any 
position in between. We induced control and Lsdl mutant clones 
in parallel. Control escort cell clones were never associated with an 
obvious robust germline tumor phenotype, although we noted one 
exception in which a single germarium with control escort cell 
clones contained six single germline cells with round fusomes (1 
out of 143 cormted). By contrast, 17% (27/155) of the germaria 
that contained Lsdl mutant escort cell clones displayed an 
expanded germline stem cell-like cell phenotype (Figure IG^. 
The relative position of l^dl mutant clones appeared to correlate 
with the appearance of a germline phenotype. The vast majority of 
germaria (96%; n = 27) that contained greater than 5 germline 
stem cell-like cells carried at least one Lsdl mutant middle escort 
cell clone. We observed one example in which a germarium with a 
mild expansion of GSC-like cells contained an Lsdl mutant 
anterior clone and posterior clone but no middle escort cell clones 
(Figure IH). Of note, most germaria that carried middle escort cell 
clones did not exhibit a GSC expansion phenotype (98/125 
germaria). WhUe their appearance was rare, germaria with only 
anteriorly or posteriorly (Figure IK) positioned escort cell clones 
did not display a robust GSC-like cell expansion phenotype. 
Similarly, loss of Lsdl in the terminal filament did not result in an 
obvious phenotype (Figure IL). 

Germaria that contained Lsdl mutant escort cell clones and 
exhibited an increased number of GSCs occasionally had an 
elongated and abnormal morphology (Figure IJ). Moreover, 22. 1 % 
of the germaria (n= 199) from Ljdl mutant females lacked marked 
escort cell clones, compared to 13.9% of control germaria (n= 166), 
and the average number of Lsdl mutant clones per germarium (4.72 
escort cell clones /germarium) was lower compared to controls (8.77 
escort cell clones /germarium), suggesting that either Lsdl mutant 
escort cell progenitors exhibited reduced proliferation during 
development or died during the course of the experiment. If 
increased death within the escort cell population accounted for all the 
observed phenotypes, one might predict that complete loss of all Lsdl 
mutant escort cell clones within a particular germarium would result 
in an increased number of GSC-like cells. However, we did not 
observe an expanded GSC phenotype in germaria that lacked clones 
from Lsdl mutant females. Together all the phenotypic data suggest 
both that escort cells require Lsdl function to limit GSC number and 
that loss of Lsdl compromises the growth and survival of escort cells, 
consistent with previous observations [24] and those noted above, 
which in turn further exacerbates the observed germ cell phenotypes. 
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Mapping Lsdl binding sites in specific cell populations 
within the germarium 

To directly define the molecular mechanisms by which Lsdl 
influences escort cell function, we elected to identify direct targets 
of Lsdl regulation in these cells. Determining whether Lsdl 
targeted potential candidate genes represented a significant 
challenge. For example, the size and complexity of the dpp 
promoter precluded our ability to assay whether Lsdl directly 
targeted this gene using a PGR based chromatin immunoprecip- 
itation (ChIP) approach. To systematically define Lsdl binding 
sites, we conducted ChIP experiments coupled with massive 
parallel sequencing (ChlP-seq). We used a number of difiFerent 
Hemagglutinin (HA) tagged transgenes, including a N-terminaUy 
tagged UASt-HA::Lsdl transgene that exhibits high expression in 
the somatic cells and a N-terminaUy tagged UASp-HAr.Lsdl 
transgene that displays relatively lower levels of expression in the 
somatic ceUs (Figure 2; SI). The UASt-HAr.LM and UASp- 
HAr.Lsdl transgenes both fuUy rescued the Lsdl^^ GSC tumor 
phenotype when driven by the c587-gal4 driver. Because Lsdl is 
expressed ubiquitously throughout the ovary [24], we sought to 
determine whether this protein bound to distinct sites in difiFerent 



cell populations within the germarium. We expressed the UASt- 
HAr.hsdl and UASp-HAr.Lsdl trangenes in cap cells and terminal 
filament cells using hh-gal4 (Figure 2 A; SI) and in the escort cells 
and early follicle cells using the c587-gal4 driver (Figure 2B; SI). 
HA-directed ChIP assays were performed on dissected ovaries and 
the immunoprecipitated chromatin was compared to input 
chromatin as a control. Within escort cells and early follicle cells, 
products of the UASp-HA::Lsdl and UASt-HA::Lsdl trangenes 
bound to 207 and 191 sites respectively (Based on the FindPeaks 
algorithm using a p-value threshold of l.OOe-3 to maximize the 
number of potential peaks; Table S1,S2), with 100 common sites 
sharing some degree of overlap (Figure 2C). Within cap cells and 
terminal filament cells, the UASp-HA::Lsdl and UASt-HA::Lsdl 
transgenic products associated with 98 and 167 genomic loci 
respectively (Table S3,S4), with 37 overlapping loci in common 
between the two datasets (Figure 2C). Comparing all four datasets 
revealed 66 common peaks between terminal filament/cap cells 
and escort cells/early follicle cells (Figure 2C,D). 232 peaks 
appeared specific for escort cells and early follicle cells and 162 
specific for cap cells and terminal filament cells (Fignre 2C,D). 
MACs analysis [27] showed similar but broader peak calls (Table 
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Figure 1 . Further characterization of the Lsd7 mutant phenotype. (A) Schematic of a Drosophila germarium. Two germline stem cells (GSCs) 
reside next to terminal filament and cap cells labeled in blue. GSCs carry a round endoplasmic reticulum-like organelle called a fusome or 
spectrosome. These fusomes become branched in appearance as GSC daughters form multicellular cysts. Lining the anterior end of the germarium 
next to the cap cells are the escort cells labeled in green. Germaria from (B) c587-gal4>mCD8::GFP and (C,D) c587-gal4>mCD8::GFP+Lsdf'^^' three day 
old females stained for Hts (red), GFP (green) and DNA (blue). RNAi knockdown of Lsdl results in the accumulation single cells with round fusomes 
and the gradual retraction of escort cell extensions. (E) c587-gal4>mCD8::GFP and (F) c587-gal4>mCD8::GFP+Lsdl"'^^' germaria from females shifted 
to 29°C for 3 days stained for activated Caspase 3 (Casp3), GFP (green) and DNA (blue). (G-L) MARCM Lsdl''" clones stained for Hts (red), GFP (green) 
and DNA (blue). Lsdl mutant escort cell clones are associated with undifferentiated germ cells with round fusomes. Scale bars= 10 jj.M. 
doi:1 0.1 371/journal.pgen.1 004200.g001 
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S5, S6, S7, S8). Lsdl enrichment peaks were spread throughout 
the Drosophila genome (Figure S2) and showed a preference for the 
promoter and 5'UTR regions of genes (Figure S3). We were 
unable to isolate a sufficient number of cells to map H3K4 
methylation across the escort cell genome. However, comparing 
our data with available data from the modENCODE project 
revealed that Lsdl binding peaks correlate with valleys of H3K4 
methylation in embryos (Figure S4), consistent with the established 
biochemical activity of Lsdl. 

Strikingly, we did not observe any enrichment for Lsdl binding 
near the dpp locus in escort cells (Figure 2E), indicating that the 
repression of BMP signal transduction by Lsd 1 must be through an 
indirect mechanism. We examined the annotation of genes near 
escort cell and early follicle cell peaks, cap cell and terminal 
filament peaks, shared peaks and from the UASt-HA::Lsdl data 
sets (Table S9, SIO, Sll, SI 2). This analysis indicated that genes 
near escort cell specific Lsdl binding peaks encode products with a 
diverse array of functions needed for development, basic cellular 
processes and reproduction (Figure S5A) [28,29]. Further analysis 
of this gene set did not reveal significant enrichment for 
components of specific pathways. MEME analysis showed an 
enrichment of ACTGGAA elements within Lsdl binding sites 
(Figure S5B). The significance of this enrichment remains unclear. 
These results suggested that the mis-expression of a functionally 
diverse set of genes likely contributes to the Lsdl mutant 
phenotypes. 

engrailed and hedgehog mis-regulation contribute to the 
Lsdl mutant phenotype 

The engraikd gene stood out as one potentially relevant target 
among the list of candidate genes, engrailed encodes a homeobox 
transcription factor that acts as a segment polarity gene during 
embryogenesis [30-32]. Prc\'i<)us results showed that engrailed 
regulates early cnarian development and that Engrailed protein 
expression is restricted to the terminal filament and cap cells in 
adult germaria [33] . Engrailed functions within these cells to help 
maintain GSCs [12]. Our ChlP-seq data revealed that Lsdl 
exhibits enriched binding to a 2 kb region of the engrailed promoter 
in the escort cells (Figure 3A). We performed RNA RT-qPCR to 
look at the transcript levels of engrailed in hidl mutants. We 
compared bam mutants to bam Lsdl double mutants because these 
samples are comparable in size and have the same basic cellular 
makeup (Figure S6). This analysis revealed that engrailed transcript 
levels increased 6 fold in the absence oi Lsdl (Figure 3B). 

Next. \\"c lcst(-d whether Engrailed protein expression expanded 
in the absence of Lsdl. In wild type germaria, cap cells and 
terminal filament cells express readily detectable levels of 
Engrailed, whereas the escort cells do not (Figure 3C) [12,33]. 
In Lsdl^''^ mutants, however, we observed ectopic Engrailed 
protein expression in a limited number of escort cells, in addition 
to the terminal filament and cap cells, in 85.7% (n = 91) of the 
germaria examined (Figure 3D). These Engrailed expressing escort 
cells tended to be positioned immediately adjacent to the normal 
niche, although occasionally we observed Engrailed expressing 
escort cells several cell positions away from the cap cells 
(Figure 3E). We cannot rule out the possibility that other cells 
also mis-expressed Engrailed, but at a level below our detection 
threshold. These data together suggest that Lsdl serves to repress 
engrailed expression within a subpopulation of escort cells. 

To test the functional relevance of ectopic Engrailed expression 
in escort cells, we assayed whether disruption of engraikd function, 
either through RNAi knockdown or loss-of-function mutations, 
modified the Lsdl mutant phenotype. Knockdown oi engrailed in an 
Lsdl RNAi background (Figure 4A) suppressed the expanded 



GSC-like cell Lsdl mutant phenotype (Figure 4B,E,F). Further- 
more, three independent engrailed mutations also suppressed the 
Lsdl RNAi-induced phenotype, so that the number of single cells 
with round fusomes decreased and cyst development within the 
germarium proceeded normally (Figure 4C-F). In all cases, 
engrailed suppression of the Lsdl RNAi ph(;notype resulted in the 
formation of morphologically normal ovarioles with maturing egg 
chambers. 

In Drosophila, Engrailed drives the expression of hedgehog (hh), 
which in turn leads to the expression of dpp in adjacent cells [34— 
3 7] . Previous analysis showed an expansion of hh expression in 
Lsdl mutant germaria [24]. To determine whether the mis- 
expression of hit in escort cells contributed to the Lsdl mutant 
phenotype, we crossed both M-specific UAS RNAi and hh mutant 
lines into a c587-gal4>UAS-Lsdl^^' background. This analysis 
revealed that loss of hh function, similar to engrailed, suppressed the 
GSC-like expansion phenotype, resulting in the formation of 
germaria that exhibited normal germ cell differentiation 
(Figure 4G-J). 

Mis-expression of engrailed in the escort cells results in a 
stem cell expansion phenotype 

To assess whether mis-expression of engrailed and hh are sufficient 
to expand the number of stem cell-like cells in the germarium, we 
expressed transgenes corresponding to both genes within escort 
cells and early follicle cells using the c587-gal4 driver. Similar to the 
phenotype exhibited by Lsdl mutants, ectopic expression of 
engrailed resulted in a stem cell-like cell expansion within 49.5% 
of germaria examined (n= 111). Many of these germline cells 
remained as single cells with round fusomes (Figure 5A). However, 
mis-expression of engrailed did not completely block cyst develop- 
ment and many ovarioles from c587-gal4>\JAS-engrailed females 
contained maturing c-gg chambers. 

In contrast to engrailed, over-expression of hh using two diflFerent 
transgenes did not obviously perturb early germ cell differentiation 
(Figure 5G,H). However, the mis-expression of these transgenes 
did result in follicle cell defects, consistent with previously 
published results [38]. These results indicated that the hh transgene 
is active in these cells but that hh over-expression in the escort cells 
and early follicle cells is not sufficient to induce an expansion of 
GSC-like cells in germaria. 

Loss of Lsdl results in expanded BMP signaling within the 
germline [24] . Based on the expansion of Engrailed expression in 
Lsdl mutants and the similarities between the Lsdl mutant and 
engrailed over-expression phenotypes, we reasoned that mis- 
expression of engrailed may also induce ectopic BMP signaling in 
the ovary. To test this, we used a Dad-lacZ enhancer trap as a 
positive transcriptional reporter oi dpp signal transduction [9,39]. 
In control ovarioles, stem cells express high levels of Dad-LacZ, 
whereas the expression of this reporter sharply decreases in 
differentiating cysts (Figure 5B). Upon engrailed mis-expression in 
the escort cells, the number of Dad-LacZ positive germline cells 
increases, likely reflecting expanded Dpp signaling (Figure 5C). 
Next, we knocked down the expression of dpp in the presence of 
the engrailed transgene and found that disruption of dpp suppressed 
the engrailed over-expression phenotype (Figure 5D-F). Together 
these results suggest that mis-expression oi engrailed in Lsdl mutants 
drives ectopic BMP signaling, resulting in the expanded GSC-like 
cell tumor phenotype. 

To test whether ectopic engraikd expression can specifically affect 
adult escort cells, we performed a temporally controlled over- 
expression experiment. c587 -gal4>\J A&-engrailed larvae were kept 
at low temperature to prevent robust expression of the engrailed 
transgene. Ovaries from adult females maintained at a low 
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temperature did not display ectopic Engrailed expression or an 
expanded undifferentiated cell phenotype (Figure 6A,A'). Howev- 
er, shifting c587-gal4>\]AS>-engrailed females to a higher temper- 
ature after eclosion resulted in ectopic engrailed expression in escort 
cells and early follicle cells, and a concomitant expansion of 
germline stem cell-like cells (Figure 6B-C). Thus, engrailed 
expression specifically in adults appears sufficient to induce ectopic 
BMP signaling in the anterior region of the germarium. 

Other Lsdl target genes 

Lsdl associates with the promoters of many genes besides engrailed, 
some of which could potentially play a role in regulating escort cell 
fimction. To begin to characterize whether Lsdl modulates the 
expression of other potential target genes, we stained c587-gal4 control 



and c587-gal4>\]ASi-L<:dl ' ' ovaries using available antibodies. Cap 
cells and escort cells exhibit a shared peak of Lsdl binding near the 
Rhol gene (Table S 1 1 , 1 2). Previous results showed that loss of Rhol 
results in escort cell defects and an expansion of GSC-Hke cells [3], 
Knocking down Lsdl levels did not appear to result in any dramatic 
change in Rho 1 expression within the escort cells (compare Figure 7 A 
and 7D). Likewise, Lsdl also exhibits enriched binding near Apcl 
(Tables S9, S 1 2), a component of the W nt signaling pathway. However 
antibody staining showed that Ape 1 protein levels did not change to an 
appreciable degree upon knock-down of I^dl (Figure 7B,E). In 
contrast, the product of a third potential Lsdl target gene. Broad, 
exhibited increased expression in c587-gal4>\]ASi-LsdlRNAi samples 
relative to controls (Figure 7C,F). However, loss of fooarfdid not appear 
to suppress the Lsdl mutant phenotype (data not shown). 
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Figure 2. Defining the targets of Lsdl in the cap cell niche and the surrounding escort cells within Drosophila germar'ia. (A,B) Germaria 
stained for HA (green), to label the transgene, Hts (red), to label fusomes, and DNA (blue) (Scale bars, 10 |ilVl). (A) hh-gal4>UASt-HA::Lsdl germarium 
showing expression in cap cells and (B) c587-gal4>UASt-HA::Lsdl germarium displaying expression in the escort cells and early follicle cells, but not in 
the cap cells. (C) Venn diagram showing the number of Lsdl binding sites in escort cells and cap cells. (D) Lsdl binding peaks on a segment of Chr2R. 
Black arrow points to a common peak in both the escort cell and cap cell populations. Orange arrows point to unique peaks seen only in the escort 
cell and early follicle cell population. (E) No Lsdl binding sites are observed around the dpp locus in escort cells. 
doi:1 0.1 371/journal.pgen.1 004200.g002 
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The raw gene functions in the developing gonad to regulate the 
morphology of somatic cells as they interact with primordial germ 
cells [40,41]. Lsdl exhibits enriched binding just 3' to the raw 
transcription termination site (Figure 7G). Antibodies were not 
available to assay whether loss of Lsdl caused a change in Raw 
expression levels but raw mutant and RNAi lines partially suppressed 
the Lsdl phenotype (Figure 7H-L). The raw''^'*''''' allele weakly 
modified the GSC4ike cell expansion phenotype exhibited by c587- 
gal4>UAS-Lsdl'^''^'^' germaria, while both raw^'^'^' and a single copy 
oi raw' ' '"^^ more strongly suppressed the c587-gal4>UAS-Lsdl^"*' 
phenotype, giving rise to a reduced number of germaria that carried 
more than 5 single cells with round fusomes. These genetic 
interactions suggest that mis-regulation of rmti expression or activity 
also contributes to the Lsdl mutant phenotype. 

Encouraged by the finding that disruption of raw suppressed the 
Lsdl mutant phenotype, we crossed a number of additional 
knockdown lines, corresponding to other putative Lsdl target 
genes, into the c587-gal4>UAS-Lsdl'^'^' background. We counted 
the total number of single germ cells with round fusomes within 
individual germaria from the resulting ovaries. This analysis 
showed that knockdown of 7 out of the 34 genes tested suppressed 
the c587-gat4> UAS-Lsdl^'^^' phenotype to the point where fewer 
than 50% of the assayed germaria contained greater than 5 single 
cells with round fusomes (Figure 8A). This group included FK506- 
binding protein 1 (FK506-bpl), Glutamirw:Jructose-6-phosphate aminotrans- 
ferase l{Gfatl), CGI 1779, ken and barbie {ken), Anaphase Promoting 
Complex submit 7 {APC7j, barren [barf] and Hepatocjte nuclear factor 4 
{Hnf4). These genes have varied functions and play roles in cell 
cycle regulation {APC7 and barr), juvenile hormone signaling 
{FK506-bpl), development of the genitalia (ken) and lipid metab- 
olism {Hnf4). Lack of a clear functional link between these 
suppressors suggests that escort cells are particularly sensitive to 
perturbations in their gene expression programs. Together these 
data show that disruption oi Lsdl results in a complex phenotype, 
marked by increased BMP signaling in the germline and 
disruption of normal escort cell function, which likely involves 
the mis-expression of several direct and potentially indirect target 
genes (Figure 8B). 

Discussion 

Using ChlP-seq, we have systematically identified Lsdl binding 
sites throughout the genomes of the cap cells/ terminal filament 
cells and the escort cells/early follicle cells of the Drosophila ovary. 
The estabhshment of ChlP-based techniques using these specific 
cell populations and this comprehensive Lsdl genome-wide 
dataset ^vill facilitate' further studies on the transcriptional 
hierarchies that regulate escort cell and cap cell function. Although 
Lsdl mutants exhibit increased dpp expression [23], Lsdl does not 
appear to directiy target the dpp locus for silencing. Rather, Lsdl 
represses the expression of the transcription factor engrailed. 
Mutations in the engrailed gene suppress the ovarian Lsdl mutant 
phenotype and mis-expression of engrailed in escort cells results in 
increased BMP signaling and an expansion of undifferentiated 
germline cells within germaria. Thus engrailed may regulate dpp 
expression in germaria through either M-dependent or indepen- 
dent mechanisms. Lsdl Eilso regulates additional genes adjacent to 
its other binding sites. Genetic analysis shows that the mis- 
regulation of these other target genes likely contributes to the Lrrfi 
mutant phenotype to varying degrees. 

Lsdl gene targets and function 

We have found Lsdl associates with a limited number of loci 
within two different cell populations. Lsdl exhibits fairly broad 



peaks of binding, ranging in size from 166-262 bp based on the 
FindPeaks algorithm. MACs analysis calls even wider peaks 
(Supplementary Tables S5, S6, S7, S8). The significance of the 
width of these peaks remains unclear but suggests that Lsd 1 either 
does not associate with single sequence specific elements at these 
sites or exhibits a certain degree of spreading upon recruitment to 
a particular locus. In Drosophila escort cells, Lsdl binds to over 100 
sites spread throughout the genome. Lsdl binds to fewer sites in 
cap cells. While some Lsdl binding sites overlap in cap cells and 
escort cells, the relatively large number of different sites suggests 
that Lsdl recruitment depends on multiple and perhaps distinct 
cell-specific co-factors. MEME analysis [42,43] reveals that many 
of the identified Lsdl binding sites contain ACTGGAA elements. 
GGAA sequences are often present in the core binding sites of 
ETS transcription factors. The Drosophila genome encodes a 
number of ETS family members, none of which have been 
characterized in the somatic cells of the germarium. Determining 
the functional relevance of these specific GGAA sites within gene 
promoters and identifying the transcription factors that bind to 
them win require additional efforts. 

For technical reasons and to enable comparisons of Lsdl 
binding between escort cells/early follicle cells and cap cells/ 
terminal filament cells, we elected to express the Lsdl HA-tagged 
transgenes in an otherwise wild-type background. We acknowl- 
edge the possibility that endogenous Lsdl may outcompete the 
HA-tagged transgenes for binding at specific sites in the escort cells 
and early follicle cells. Therefore sites identified in this study may 
be an underrepresentation of the total number of sites bound by 
endogenous Lsdl. Repeating the ChlP-seq analysis using material 
from rescued Ljdl^^^ females that express the HA-tagged Lsdl 
transgene in escort cells and early follicle cells represents important 
work for the future. 

We found that Lsdl mutant samples exhibit a 6-fold increase in 
engrailed transcript levels relative to controls. Curiously, ectopic 
Engrailed protein expression was only observed in a small number 
of escort cells. Perhaps additional post-transcription mechanisms 
regulate the translation of Engrailed, and potentially other proteins, 
inside and outside of the cap cell niche. Such mechanisms would 
allow these cells to fine-tune their signaling output more than what 
could be achieved through transcriptional based mechanisms alone. 
Results presented here also suggest that escort cells are not uniform 
in nature and perhaps have specific functions or capabifities 
depending on their lineage and where they reside within the 
germarium. MARCM analysis shows that the loss of Lsdl in some 
but not all escort cells results in a marked expansion of GSC-like 
cells within the germarium. Previous studies have also suggested that 
specific escort cells have distinct roles in supporting GSCs [12]. 

Technical considerations aside, the severity of the Lsdl nuU 
phenotype compared to the engrailed over-expression phenotype, 
both in terms of the penetrance and extent of the GSC expansion 
phenotype and the accompanying germline and somatic cell 
death, suggests that engrailed is not the only biologically relevant 
target of Lsdl regulation in the escort cells. Based on expression 
analysis (Figure 7), Lsdl regulates the expression of some but not 
all genes adjacent to its binding sites. Our genetic analysis suggests 
that mis-regulation of the putative target raw (Figure 7) and several 
additional genes (Figure 8) also contribute to the Lsdl mutant 
phenotype. Characterizing the transcriptional profile of escort cells 
from wild-type and L^dl mutant samples, which will have to await 
for improvements in current cell isolation and RNA profiling 
techniques, will help to further resolve which genes are direct and 
indirect targets of Lsdl regulation. Such approaches may also 
reveal additional genes that participate in niche formation and 
function. Of note, the observation that functionally diverse genes 
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Figure 3. Lsdl silences engrailed'in the escort cells. (A) A Lsdl binding peak is observed in the engrailed promoter in escort cells. (B) RT-qPCR of 
whole ovary extracts from bam'^^^ and bam'^^^ LsdT"^ double mutants. The bam^^^ Lsdl^'^ double mutant samples exhibit a relative 6-fold increase of 
engrailed transcripts when com pared to bom""*'^ mutants. (C and C) w'"" and (D-E') Lsdi^^ homozygous germaria stained for Engrailed (En) (green) 
and VASA (red) and DNA (blue). Arrows point to Engrailed positive cells. In w'"^ control germaria, Engrailed is expressed in the terminal filament and 
cap cells, whereas Lsdi mutants display engrailed expression in the anterior escort cells. Scale bars = 10 |iM. 
doi:1 0.1 371 /journal.pgen.1 004200.g003 



can suppress the Lsdl mutant phenotype suggests that escort cells 
are acutely sensitive to changes in their gene expression programs. 
While our data support a model that loss of Lsdl initially results in 
mis-regulation of engrailed and other genes that, in turn, drive GSC 
expansion, it is clear that many of the escort cells that experience 
reduced Lsdl function retract their cellular extensions and 
undergo cell death. This loss of escort cells further exacerbates 



the GSC expansion phenotype. Given the phenotypic complexity 
described here and elsewhere [3], care should be taken when 
analyzing gene function within the escort cell population. 

Niche plasticity 

How niches maintain stem cells and adjust their signaling 
output to ensure tissue homeostasis remains a fundamental 
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Figure 4. engrailed mutations suppress tfie Lsd7 RNAi pKienotype. Ovarioles from (A) c587-gat4>UAS-Lsdf"'^' (B) c587-gal4>UAS-L5dl""^'; 
UAS en™*'"V+ (C) c587-gal4>UAS-Lsdf"^'; en7+ and (D) c587-gal4>UAS-Lsdf^^'; enV+ stained for Hts (green), VASA (red) and DNA (blue). (E) Graph 
showing the percentage of germaria that contain the indicated number of single germ cells with round fusomes for each genotype. (F) Graph 
showing the percentage of ovarioles that have a given number of developing egg chambers for each genotype. c587-gal4> UAS-Lsdl'^'^'^' (n = 171); 
c587-gal4>UAS-Lsdf"^'; en7+ (n = 90); c587-gal4>UAS-Lsdf'^^'; en""/+ (n = 95); c587-gal4>UAS-Lsdf'^^'; enV+ (n = 115); c587-gal4>UAS-Lsdf"''''; 
6^7"™'-'/+ (n = 108); c587-gal4>UAS-Lsdf"^'; en''™'-7+ (n = 96). (G) c587-gat4>UAS-Lsdr^''''; UAS (H) c587-gat4>UAS-Lsdf^'^'; hh'"^/+ and (1) 

c587-gal4>UAS-Lsdf"'^'; hh^/+ stained for Hts (green), VASA (red) and DNA (blue). (J) Graph showing the percentage of germaria that contain the 
indicated number of single germ cells with round fusomes for each genotype. c587-gal4>UAS-Lsdl™^' (n = 307); c587-gal4>UAS-Lsdf"^'; UAS h/i™'^'/ 
+ (n = 222) c587-gal4>UAS-Lsdl™'^'; hh^^/+ (n = 88) and c587-gal4>UAS-Lsdf"^'; hh^/+ (n = 184). Scale bars =10 |ilVI. 
doi:10.1371/journal.pgen.1004200.g004 



question in stem cell biology. Elegant work has shown that 
terminal filament cells, cap cells and escort cells help to support the 
self-renewal of two to three germline stem cells at the tip of 
Drosophila germaria [5,44]. The predominant signal emanating 
from the anterior tip of the germarium is Dpp, which acts locally 
to induce a canonical signal transduction cascade in GSCs, which 
in turn represses their differentiation [5,9] . Several expression and 



genetic studies strongly suggest that terminal filament and cap 
cells, and perhaps the most anterior escort cells, are the primary 
source Dpp ligand [5,9]. More recent work has suggested that 
Engrailed expression in cap cells non-autonomously promotes dpp 
expression in escort cells through a hedgehog dependent mechanism 
[12]. Loss oi Lsdl results in ectopic expansion of expression in 
escort cells [24] and data shown here (Figure 4) reveals that 
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Figure 5. Ectopic expression of engrailed in the escort cells results in an expansion of stem cell-like cells in the germline. (A) c587- 
gal4>UAS-en::GFP germaria stained for Hts (green), VASA (red) and DNA (blue). (B) Control and (C) c587-gal4>UAS-en::GFP germaria stained for Dad- 
LacZ (Dad) (green), VASA (red) and DNA (blue). In control germaria, Dad-LacZ is expressed at high levels in the germline stem cells and its expression 
is reduced in the early differentiating cells. Germaria overexpressing engrailed display an expansion of Dad-LacZ positive cells. (D) Germaria from 
c587-gal4>UAS-en::GFP, UAS-dpp"'^'"'^ and (E) c587-gal4>UAS-en::GFP, UAS-dpp""'^" ovaries stained for Hts (green), VASA (red) and DNA (blue). 
Reducing the levels of dpp in the escort cells suppressed the engrailed overexpression phenotype. (F) Graph showing the percentage of germaria that 
contained a given number of cells with round fusomes for each genotype. (G) c587-gal4>UAS-hh(lll) and (H) c587-gal4>UAS-hh(ll) stained for Hts 
(green), VASA (red) and DNA (blue). Scale bars = 10 nM. 
doi:1 0.1 371 /journal.pgen.1 004200.g005 



disruption of hh partially suppres.ses the Lsdl mutant phenotype. 
However, consistent with previous results [38], over-expression of 
hh in escort cells does not result in an Lsdl-like mutant tumor 
phenotype (Figure 5G,H), demonstrating that hh is not sufficient to 
induce ectopic BMP signaling in the germarium. Given these 
observations, ectopic mgraikd expression in escort cells likely 
targets additional genes besides hh to induce ectopic BMP signaling 
and promote the expansion of undifferentiated germ cells. 

The finding that loss of Lsdl or mis-expression of engrailed in 
adult escort cells leads to expanded Dpp signal transduction within 
germ cells throughout the anterior portion of the germarium 
indicates that subpopulations of escort cells are capable, and 
perhaps even poised, to express dpp under certain conditions. Such 
plasticity might allow the niche to expand and contract in response 
to various stimuli and environmental cues. Indeed, previous 
studies have shown that Jak/Stat and insulin signaling can 
influence the number of GSCs in the ovary [45-49] . Moreover, 
ongoing dynamic regulation of signaling may be a regular feature 
of niches under resting homeostatic conditions. The observation 
that long-term knock-down of dpp in escort cells results in a 
reduced number of GSCs at the tip of the germarium, but not 
their complete elimination, is consistent with the notion that escort 
cells contribute to the maintenance of GSCs in some manner [12]. 
Further work, with single-cell spatial and small-scale temporal 
resolution, will be needed to help clarify what cells express niche 
signals and when. 

Inappropriate and extensive expansion of niches would be 
predicted to upset tissue homeostasis and perhaps even result in 
pathological conditions. Therefore robust but flexible mechanisms 
that depend on chromatin factors such as Lsd 1 may be in place to 



precisely control the expansion and contraction of in vivo stem cell 
niches. The continued study o{ Dro.sophila cap cells and escort cells 
will provide further insights into how chromatin programming 
regulates niche plasticity. 

Materials and Methods 

Drosophila stocks 

Drosophila stocks were maintained at room temperature on 
standard cornmeal-agar medium unless specified otherwise. The 
foUowing fly strains were used in this study: w'"^ was used as a 
control; Lsd 1'^'^ was provided by N. Dyson (Massachusetts General 
Hospital Cancer Center, Charlestown, MA); hh-gal4 and UAS-M 
lines [50] were provided by J. Jiang (University of Texas 
Southwestern, Dallas, TX); c587-gal4 and Dad-Lac^wexe. provided 
by A. Spradling (Carnegie Institution for Science, Baltimore, MD); 
the UAS-engrailed::GFP transgenic line was provided by Florence 
Maschat (Institute of Human Genetics, France); the raw''^^'"*^ and 
raiu^^''^'' alleles were provided by Jennifer Mierisch (Loyola 
University of Chicago); er^ , m , en^', UAS-mCD8::GFP, UAS- 
dpp'^'^'-' (BL#-31172), UAS-dpp^^''^ (BL#-31530), UAS-dpp^-*''^^ 
(BL#-31531) and UAS-raw^-*' {BL#-31393), UAS^en^"^' ' (BL#- 
33715), (BL#-26752), K45'-M^^'^^' (BL#-31042), M^'^ 

(BL#-1749), hh^ (BL#-3376), hroad''''^ (BL#-29971), Aroo/ (BL#- 
29972), par-l''^''""*'" (BL#- 32410), NaPi-l"'^'""'""' (BL#- 
34003), CGniS^""'^*'^'^ (BL#- 27079), CG12054"'^"^'^'^ (BL#- 
50910), Atgl""™'-'^^" (BL#- 44034), muS^^'^""'-'" (BL#- 35044), 
CGI 2 128"^^^""'""' (BL#- 33997), Nek2"^'^"^"^" (BL#- 28600), 
GstSl"^^"-'^'"''^ (BL#- 28885), RhoGEF2"^'^"" (BL#- 34643), 
1^^1.^0164!, ^gL#- 37506), UGP"'^^"'"'"' (BL#- 50902), 
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Figure 6. engrailed mis-expression in adult escort cells results in delayed germ cell differentiation. c587-gal4>UAS-en::GFP feinales were 
kept at (A) 18°C during both development and in adulthood or (B) moved to 29'C after eclosion. (A,B) The germaria were stained for alpha Spectrin 
(Spec) (green), VASA (red) and DNA (blue) and (A',B') Engrailed (En) (grayscale). (C) Graph showing the percentage of germaria with the indicated 
number of single germ cells with round fusomes when subjected to different temperatures. c587-gal4>UAS-en::GFP 18°C>18"C (n = 100), c587- 
gal4>UAS-en::GFP 18 C>29°C (n = 80). Scale bars = 10 ^iM. 
doi:1 0.1 371 /journal.pgen.1 004200.g006 
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(BL#- 34068) and //n/-^^"^'^" (BL#- 29375) lines 
were obtained from the Bloomington Stock Center. UAS-Lsdl^^^ 
was obtained from the National Institute of Genetics, Japan. 

MARCM clones 

Lsdl mutant MARCM clones were generated by crossing 
Udl"^^ FRT 2A to yw'^^ tub-gal4 UAS-GFP;;tub-gal80 FRT 2A/ 
TM6B (gift from Ben Ohlstein). The resulting larvae were heat- 
shocked twice a day at 37°C on days 5-7 after the cross was set. 
The resulting adult females were dissected and stained 7 days after 
they eclosed. 

Cloning of tagged transgenes 

The HA tagged transgenes of Lsdl were created using Gateway 
Cloning (Invitrogen). The open reading frame (ORF) of Lsdl 
was cloned into modified pTHW and pPHW destination 
vectors (http:// emb.carnegiescience.edu/labs/murphy/ Gateway% 



20vectors.html) that contained (pC31 attB sites [51-53]. These 
constructs were injected into flies and transformed using (pC31 
integrase into die 5 ID landing site on the 2"'' chromosome. 

Immunostalning 

Adult ovaries were dissected in Grace's medium and fixed in 
4% (vol/vol) formaldehyde for 10 mill. The ovaries were washed 
with PBT (IX PBS, 0.5% BSA, and 0.3% Triton-X 100) and 
stained with primary antibody overnight at 4°C. The ovaries 
were washed and incubated in secondary antibody at room 
temperature for 5 hrs. Ovaries were then washed again and 
mounted in Vectasliield containing DAPI (Vector Laboratories). 
The following primary antibodies were used: mouse anti-Hts 
(IBl) (1:20), mouse anti-Engrailed (4D9) (1:2), mouse anti-Broad- 
core (25E9.D7) (1:10), mouse anti-Rhol (P1D9) (1:50) and rat 
anti-VASA (1:20) (Developmental Studies Hybridoma Bank), 
mouse anti-(3-galactosidase (1:200) (Promega), rabbit anti-APCl 
[54] (l:1000)(gift of E. Wieschaus), Rabbit anti-a-Spectrin 
[55](l:1000)(gift from Ron Dubreuil), rabbit anti-GFP (1:1000) 
(Molecular Probes), rat anti-HA 3F10 (Roche) and rabbit anti- 
cleaved Caspase-3 (1:250) (Cell Signaling Technology). Fluores- 
cence-conjugated secondary antibodies Jackson Laboratories) 
were used at a dilution of 1:200. 

RNA isolation and RT-PCR 

RNA was isolated from bamA86 and Lsdl^''^' bamA86 mutant 
ovaries using TRIzol (Invitrogen). The RNA was treated with 



PLOS Genetics | www.plosgenetics.org 



10 



IVIarch2014 | Volume 10 | Issue 3 | e1004200 



Lsdl Regulates Escort Cell Function 



c587-gal4 



c587-gal4>Lsd1 



RNAi 



A 

— s 


. ■ u 

pec DNA 


— Rhol 


D 

— Rhol 


EjI 




APC1 




Hts Vdsa DNA 


c 

f ■ 

— s 


pec DNA 


Br 




• • •• 



G 

3 

c587>Lsd1 

c587>lnput 
hh>Lsd1 



100 



CG9314 



Hnf4 



IIIHI 



H c587>UAS-Lsdl'"''' 
_ H»\/?fsa DNA 


1 cS87>UAS-Lsd1'""";raw"V+ 
* . V/- 

_ HtsVasaDNA 


J c587>UAS-LsdT<'^';rav/'"l+ 


K c587>U^S-/.sc/7'""";raiv«™'/+ 




_ HtsV/asaDNA 



cS87-gal4>UAS-Lsd1""''' Sm^le 




Figure 7. raiy represents another potential functional Lsdl target gene. (A-C) c587-gal4 and (D-F) c587-gal4>UAS-Lsdl'^'^^' germaria stained 
for Rhol (A,D), APC1 (B,E), and Broad (Br) (C,F). Alpha-Spectrin (Spec) or Hts (green), Vasa (red) and DNA (blue) staining of the same germaria are 
shown in color. (G) Screen shot of raw locus showing a peak of Lsdl binding in the c587-gal4>UASt-HA-Lsdl sample relative to the input DNA. (H) 



c587-gaH>UAS-Lsdf"'" , (I) c587-gal4>UAS-Lsdl™^'; UAS raw'^y+ (J) c587-gal4>UAS-Lsdf 



■V+ and (K) c587-gal4>UAS-Lsdf 



stained for Hts (green), VASA (red) and DNA (blue). (L) Graph showing the percentage of germaria that contain the indicated number of single germ 
cells with round fusomes for each genotype. c587-gal4>Lsdf'''^' (n = 351 ); c587-gal4>UAS-Lsdf'^^'; UAS rav\/"'"^'/+ (n = 1 55); c587-gal4>UAS-Lsdf'^^'; 



raw"'^'"/+ (n = 224); and c587-gai4>UAS-Lsdf'^'^'; 
doi:1 0.1 371 /journal.pgen.1 004200.g007 



raw""V+ (n = 229). Scale bars = 1 0 ixM. 



DNase and subjected to RT-qPCR reaction using the Superscript 
III First-Strand Synthesis SuperMix (Invitrogen). The primers 
used to amplify engrailed mRNA are as follows: 

en^raiW forward: 5' - GCCCGCCTGGGTGTACTG 
engrailed reverse: 5' - CGCTTCTCGTCGTTGGTCTTG 

Chromatin immunoprecipitation and sequencing (ChlP- 
seq) 

We used 1000 pairs of ovaries for each ChlP-seq reaction. 
Every 200 pairs of ovaries were dissected, fixed, washed and 
frozen immediately at — 80°C. The entire protocol was done at 



4°C unless otherwise indicated. The ovaries were dissected and 
fixed in 1 ml of 1% formaldehyde at room temperature for 
10 mins. The crosslinking was stopped by the addition of 100 ml 
1.25M glycine solution. The ovaries were washed three times 
with IX cold PBS buffer and then sonicated in 500 |xl ChIP 
Sonication Buffer (l%Triton X-100, 0.1% Deoxycholate, 50 mM 
Tris 8.1, 150 mM NaCl, 5 mM EDTA) on ice to achieve a final 
DNA length of 100 to 600 base pairs. The sample was centrifuged 
at maximum speed at 4°C to remove cell debris. The supernatant 
was transferred to a new tube and the sonicated sample was then 
blocked by adding Protein G agarose beads and incubating at 
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Figure 8. Further genetic analysis shows that reduced levels of several genes suppress the Lsdl mutant phenotype. (A) The graph 
shows the percentage of germaria that contain the Indicated number of single germ cells with round fusomes for each genotype. Control 1 
represents c587-gal4>UAS-Lsdf"^'; +/+ and the rest represent c587-gal4>UAS-Lsdr"^'; gene specific RNAi/+. The resulting females were shifted to 
29°C for seven days. The number (n) of germaria counted for each sample Is Indicated above the bar. The yellow dotted line marks the 50% threshold. 
All the Indicated lines were crossed on their own to the c587-gal4 driver. None of the resulting ovaries displayed phenotypes except for 
CGI 3192"'^^°'^^ (egg chamber defects) and par-l""^""""" (follicle cell defects). (B) Model of Lsdl function In escort cells. 
dol:10.1371/journal.pgen.1004200.g008 



4°C for one hour. The beads were removed. 1 % of the sample 
was kept aside as INPUT and to the remaining sample 3 ug 
rabbit-HA antibody (Abeam) was added and incubated overnight 
at 4°C. 

The next day protein agarose G beads were added and 
incubated for 3 hours at 4°C. The beads was then washed well 
with Chip Sonication Buffer (two times), High Salt Wash Buffer 
(1% Triton X-1 00, 0.1% Deoxycholate, 50 mM Tris 8.1, 500 mM 
NaCl, 5 mM EDTA) (three times), LiCl Immune Complex Wash 
Buffer (two times) and TE buffer. The protein bound to the beads 
was eluted using 500 \A Elution Buffer (1% SDS, O.IM NaHCOg). 
The elution buffer was added to the INPUT samples and they 



were treated the same as the IP samples from this point. 20 [ll of 
5M NaCl was added to 500 |al of elution buffer and incubated at 
65°C overnight. 

The third day, the sample was treated with RNase A, Proteinase 
K and the DNA isolated using Qiagen PCR Purification kit. 
Subsequent library construction and sequencing of the input and 
immunoprecipitated DNA were conducted by the UT Southwest- 
ern McDermott Sequencing Center. 

Bioinformatic analysis of ChlP-Seq data 

The primary sequencing data was mapped to the fly reference 
genome dm3 using BioScope (1.2.1). During the alignment, three 
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filter steps were applied to remove low quality, ambiguous and 
redundant reads. HA-Lsdl binding regions were identified as 
genomic regions with a significant read enrichment and binding 
peak profile in the ChIP reads over the input reads using the 
FindPeaks module in the Homer software tool [56] with 10% 
false discovery rate (FDR). ChIP enrichment at important 
genome features such as specific chromosomes, promoters, 
downstream, exonic, intronic and distal intergenic regions was 
statistically analyzed with the Cis-regulatory Element Armotation 
System (CEAS) [57]. De novo motif discovery analysis for HA- 
Lsdl binding regions was performed with the Multiple EM for 
motif elicitation (MEME) software tool [42,43]. High quality 
motifs were aligned against transcription factor motifs retrieved 
from JASPAR [58] and TRANSFAC [59] using tire TOMTOM 
software tool [42] to identify known transcription factor motifs 
that match the MEME predicted motifs. Potential protein-coding 
target genes associated with the identified HA-Lsdl binding 
regions wer(; id(;iitifi(;d based on the distance of their transcrip- 
tion start sites (TSSs) according to their RefSeq annotation in the 
dm3 assembly to binding peak summits. Genes with TSSs within 
5 kb or nearest to an FL^-Lsd peak summit were called as target 
genes. ChlP-Seq data has been deposited with NCBI GEO 
(http://www.ncbi.nlm.nih.gov/geo) under the accession code 
GSE54376. 

Supporting Information 

Figure SI Expression patterns of different Lsdl transgenes when 
driven by cap cell or escort cell specific drivers. AU the images were 
captured using the same confocal settings so the relative expression 
levels could be compared between the samples. (A,B,C,D) 
Germaria stained for HA (green), Hts (red) and DNA (blue). 
(A',B',C',D') HA staining alone. (A) c587-gd4>UASt-HA::LM 
and (B) c587-gal4>UASp-HA::Lsdl germaria show expression in 
escort cells and early follicle cells. (C) hh-gal4> UASt-HA::Lsdl and 
(D) hh-gal4>UASp-HA::Lsdl germaria display expression of 
HA::Lsdl in the cap cells and the occasional terminal filament 
cell. (Scale bars, 10 \xM). 
(TIF) 

Figure S2 The chromosomal distribution of HA::Lsdl binding 
sites identified using FindPeaks in c587 -gal4> UASt-HA::Lsdl , c587- 
gal4>UASp-HA::Lsdl, hh-gal4>UASt-HA::Lsdl and hh-gal4>UASp- 
HAr.Lsdl germaria. The blue bars represent the percent of the 
genome comprised by each chromosome feature while the red 
bars indicate percentage distribution of Lsdl binding sites across 
each chromosome feature. This analysis reveals that, while 
hh>pTHW displays a modest enrichment for binding on the 
right arm of chromosome 3, in general, Lsdl binding appears 
evenly distributed across the genome. 
(TIF) 

Figure S3 The distribution of HA::Lsdl binding sites relative to 
gene features in c587-gal4>UASt-HA::Lsdl , c587-gal4>UASp- 
HA::Lsdl, hh-gal4>UASt-HA::Lsdl and hh-gal4>UASp-HA::Lsdl 
samples. 
(TIF) 

Figure S4 Read depth plots for input DNA, H3K4me3 and 
H3K4mel ChIP DNA from die modENCODE project within +/ 

— 3 kb of superimpo.sition of all Lsdl binding sites (defined as 0). 
Valleys of H3K4me levels exist in regions corresponding to escort 
cell Lsdl binding sites. 
(TIF) 

Figure S5 (A) Gene ontology hierarchies of the next gene 
adjacent to escort cell specific Lsdl binding sites (no distance 



cutoff) or those genes with transcriptional start sites within 2.5 kb 
of Lsdl binding sites based on the UASt-HA::Lsdl data sets 
[28,29]. Yellow to Orange represents less significant to more 

significant terms. The size of nodes corresponds to the number of 
genes in the query set that belong to the category. (B) A motif 
enriched in Lsdl binding sites detected by MEME analysis. 
(TIF) 

Figure S6 (A) barn^'^'' and (B) barn^'^'' Lsdl^^ double mutant 
germaria stained for Hts (green), Vasa (red) and DNA (blue). Both 
the single and double mutant germaria are roughly the same size 
and comprised of similar cell types. Scale bars= 10 nM. 
(TIF) 

Table SI FindPeaks output for c587-gal4> UASp-HA::Lsdl ChlP- 

seq. 

PCLSX) 

Table S2 FindPeaks output for c587-gal4> UASt-HA::Lsdl ChlP- 
seq. 

PCLSX) 

Table S3 FindPeaks output for hh-gal4>UASp-HA::Lsdl ChlP- 
seq. 

PCLSX) 

Table S4 FindPeaks output for hh-gal4>UASt-HA::Lsdl ChlP- 
seq. 

PCLSX) 

Table S5 MACs output for c587-gal4>UASp-HA::Lsdl ChlP- 
seq. 

PCLS) 

Table S6 MACs output for c587-gal4> UASt-HA::Lsdl ChlP-seq. 

(XLS) 

Table S7 MACs output for hh-gal4>UASp-HA::Lsdl ChlP-seq. 
PCLS) 

Table S8 MACs output for hh-gal4> UASt-HA::Lsdl ChlP-seq. 

(XLS) 

Table S9 Genes near escort cell and early follicle cell peaks. 
pCLSX) 

Table SIO Genes near cap cells and terminal filament cells. 

(XLSX) 

Table Sll Genes near shared peaks. 
pCLSX) 

Table S12 Genes within 5 kb of UASt-HA::Lsdl binding peaks. 

(XLSX) 
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